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The need for a local source of UHE CR nuclei 
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Recent results of the Pierre Auger (Auger) fluorescence detectors indicate an increasingly heavy 
composition of ultra-high energy (UHE) cosmic rays (CRs) . Assuming that this trend continues up 
to the highest energies observed by the Auger surface detectors we derive the constraints this places 
on the local source distribution of UHE CR nuclei. Utilizing an analytic description of UHE CR 
propagation we derive the expected spectra and composition for a wide range of source emission 
spectra. We find that sources of intermediate-to-heavy nuclei are consistent with the observed 
spectra and composition data above the ankle. This consistency requires the presence of nearby 
sources within 60 Mpc and 80 Mpc for silicon and iron only sources, respectively. The necessity of 
these local sources becomes even more compelling in the presence nano-Gauss local extragalactic 
magnetic fields. 

PACS numbers: 13.85.Tp, 98.70.Sa 



I. INTRODUCTION 



Ultra-high energy CRs with energies above 10"'^^ eV 
arrive at Earth with a frequency of less than one event 
per square-kilometer-year (i.e. with an energy flux of 
30 eV cm~^ s""'^) in tt steradian. The recently completed 
Pierre Auger observatory [1 with a surface detector cov- 
ering an area of about 3000 kni^ is thus able to detect 
up to several hundred UHE CR events per year. How- 
ever, due to the steeply falling CR spectrum the arrival 
frequency drops by about two orders of magnitude as we 
go up in energy by one decade, leaving only a few events 
per year detectable at energies around 10^° eV. Hence, 
the statistical uncertainty associated with the upper end 
of the spectrum is still large, limiting our knowledge of 
UHE CR composition and origin. 

Before their arrival, UHE CRs must propagate across 
the astronomical distance between their source and 
Earth. The relevant interactions of UHE CR nuclei dur- 
ing propagation are Bethe-Heitler pair production and 
photo-disintegration in collisions with photons of the cos- 
mic background radiation. The cross-sections of both 
these processes rise quickly above the threshold values of 
about 1 MeV and 10 MeV, respectively, in the nuclei's 
rest-frame. At even higher energies above 150 MeV pion 
production turns on and becomes the dominant energy 
loss process. However, for the UHE CR cutoff energies 
and composition we consider, which are motivated by 
the most recent Auger results, this process never plays a 
dominant role and may be safely neglected. 

Provided the propagation time from their sources 
to Earth is greater than their energy loss time, UHE 
CRs invariably undergo these energy loss interactions. 
The break-up of nuclei via photo-disintegration produces 
lower mass nuclei with the same Lorentz factor. For 
heavy nuclei, the most dominant transitions are one- 



nucleon and two-nucleon losses. Secondary heavy nuclei 
remain close to the photo-disintegration resonance and 
quickly disintegrate further to lighter nuclei. Hence, on 
resonance, the initial mass composition of the sources 
is quickly shifted to lower atomic number values and in 
general shows a strong dependence on CR energy. 

The arriving flux from an ensemble of UHE CR sources 
can be expected to contain suppression features at the 
high energies at which the photo-disintegration processes 
turn on and the nuclei particle's attenuation length de- 
creases. Most prominently, for the case of a proton- 
dominated spectrum, the flux is expected to be sup- 
pressed by the Greisen-Zatsepin-Kuz'min (GZK) cut- 
off [H [3] due to resonant pion photo-production inter- 
actions with the cosmic microwave background (CMB). 
Intriguingly, a suppression of the CR spectrum at about 
5 X 10^^ eV has been observed at a statistically signif- 
icant level [H [S] . However, since a similar feature may 
also appear in the spectrum from nuclei primaries, the 
observation of such a feature provides little clue as to 
the underlying composition. Such a suppression feature 
becomes a cutoff in the arriving flux if the attenuation 
length drops below the distance to the nearest UHE CR 
source. Thus, the shape of the suppression/cutoff feature 
does contain information about the source distribution, 
as has been investigated already for the case of UHE CR 
protons HHS]. 

On their arrival at Earth, UHE CRs interact with 
molecules in the atmosphere and deposit their energy in 
the form of extensive air showers. The characteristics 
of these showers along the shower depth X (in g/cm^) 
contain vital UHE CR composition information. On av- 
erage, proton-induced showers reach their maximum de- 
velopment, (^max), deeper in the atmosphere than do 
showers of the same energy generated by heavier nuclei. 
Accompanying this effect, the shower to shower fluctu- 
ation of ATmax about the mean, RMS(Amax), is larger 



for proton-induced showers than for iron-induced show- 
ers of the same energy. As a result, measurements of both 
(-'^max) and RMS(Xinax) can be used to infer the average 
chemical composition of the UHE CRs as a function of 
energy. 

Auger has now released their first measurement re- 
sults of RMS(Xmax) 011], along with those of (X^ax)- 
These results seem to imply that the UHE CR spec- 
trum contains a large fraction of intermediate or heavy 
mass nuclei, becoming increasingly heavy at high en- 
ergies (~ 10^^'^ eV). Furthermore, the small values of 
RMS(A'niax) measured by Auger also imply that UHE 
CRs are composed of species with a relatively narrow 
distribution of charge at the highest measured energies, 
containing little or no protons or light nuclei. In this way, 
the new RMS(A'inax) measurements not only confirm and 
reinforce the conclusions drawn from their earlier aver- 
age depth of shower maximum measurements, but also 
provide complementary information that enables one to 
constrain the distribution of the various chemical species 
present within the UHE CR spectrum. Curiously, re- 
cent corresponding measurements by the Telescope Ar- 
ray UHE CR detector do not appear to agree with the 
Auger results P^ , and instead seem to be consistent with 
a light composition as was indicated by its predecessor 
HiRes [TT] . Though an understanding of this conflict is 
crucial for future progress in the field, we here chose to 
adopt only the Auger results. 

Interestingly, the intermediate-heavy nuclei scenario 
emerging from these measurements may find consistency 
with the presently observed mild departure from isotropy 
and absence of point sources observed at UHE [13] . This 
could occur through the introduction of ~ 10° deflections 
in the turbulent (/iG) Galactic magnetic field structure 
[Hj or in 0.1 uG turbulent extragalactic field structure 
[15j . from an anisotropic distribution of nearby sources. 
Such large deflections would limit the prospects for fu- 
ture UHE CR astronomy for all but the brightest and 
closest sources. The correlation of UHE CR with nearby 
objects, however, still awaits to be resolved within an 
intermediate-heavy nuclei framework, and deserves fur- 
ther consideration. 

Following our previous investigations into the composi- 
tion of UHE CRs lE\, an intermediate- heavy (silicon-to- 
iron) type composition was found to be motivated by the 
present complete Auger data set, with a hard (a < 2) in- 
jection spectral index [45]. Furthermore, this agreement 
received only very mild improvements by the additional 
consideration of an admixture of species. Building fur- 
ther on these results, we here focus on the local UHE CR 
source distribution inferred to exist if the observed trend 
in the composition continues up to highest energies ob- 
served by the ground array (~ 10^°'^ eV). By considering 
simple silicon and iron source composition scenarios, we 
investigate how far away these sources can afford to be 
without detrimentally damaging the agreement with the 
Auger data. To aid this investigation we utilise an ana- 



lytic description of UHE CR nuclei propagation already 
developed by the authors [17]. 



II. UHE CR NUCLEI FLUXES FROM A 
UNIFORM DISTRIBUTION OF SOURCES 



With no prior knowledge about the UHE CR source 
population, a "universal" (homogeneous and isotropic) 
distribution is generally assumed. We adopt this assump- 
tion here as a means of investigating signatures of a de- 
parture from it. In order to quantify the effect of a differ- 
ent source distribution in this paper, we separate out the 
fluxes produced from source regions with shells of radii 0- 
3 Mpc, 3-9 Mpc, 9-27 Mpc, 27-81 Mpc, and 81-243 Mpc 
surrounding the Earth. In this way, the results obtained 
may be used to encapsulate the effects introduced by a 
non- "universal" local void of UHE CR sources. 

We show in the upper-panel of Fig. [l] a breakdown of 
the total arriving flux from an ensemble of sources as- 
suming they have a homogeneous distribution and emit 
a purely iron type UHE CR nuclei composition. These 
results have been obtained using both a Monte Carlo 
description of UHE CR nuclei propagation, whose de- 
tails are described in Ref. [TH], as well as an analytic 
description whose details are described below. An in- 
jection spectrum of a = 1.8 and a maximum energy of 
Epcmax — 10^^ eV have here been assumed [35], We use 
the iron equivalent value of the cutoff, -E'pe.max, as the 
reference value since this allows a rigidity independent 
comparison of the cutoffs values for the type species we 
consider. The injection spectrum values adopted in Fig.[T] 
were motivated by previous investigations into UHE CR 
nuclei sources [16] . From this flgure, it is clearly seen 
that the arriving flux at the highest UHE CR energies 
{E > 10^° eV) is dominated by the local source popula- 
tion. 

In the middle and bottom panels of Fig.jljthe composi- 
tion related information of the arriving flux, (Xmax) and 
RMS(Xiuax)i are shown. For these results, the hadronic 
interaction model QGSJET 11 [^0] has been adopted. 
The broad red and blue lines depicting the "proton" and 
"iron" values in these plots indicate the spread of pre- 
dicted values from a range of hadronic interaction mod- 
els [2U^E5] . The departure of the black line from the blue 
region in both these panels thus demonstrates the degree 
of photo-disintegration incurred by the injected iron nu- 
clei en route in the extragalactic radiation fields. 

To simplify the picture the effects of the presence of 
extragalactic magnetic fields on UHE CR propagation 
are neglected in this section. In section [TV] these effects 
are incorporated to see how the field-free results in this 
section are altered. 

The shape of the cutoff feature in the UHE CR spec- 
trum carries within it valuable information about the lo- 
cal source distribution. A simple analytic description of 
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FIG. 1: A comparison of the Monte Carlo and analytic flux 
results following the emission of Fe nuclei from sources with 
i?Fo,max ~ 10^^ cV and a — 1.8, to recent Auger measure- 
ments [9l 1121 119) . Top panel: A breakdown of the arriv- 
ing flux from shells of UHE CR Fe sources. Middle and 
bottom panels: The arriving flux's corresponding (Xmax) 
and RMS(Xinax) values, which describe the composition. The 
hadronic model QGSJET 11 [50] has been adopted for these 
comparative calculations. The red and blue bands show the 
range of predictions of the {X^ax) and RMS(Xniax) values 
for protons and iron, for various hadronic interaction models 
|20H23] . In all three panels the arrows labeled asya depict the 
size of the Auger systematic errors. 



[T7] , allows the different cutoffs in the fluxes from the var- 
ious shells to be easily interpreted. To first order, photo- 
disintegration of heavy nuclei can be approximated by 
one-nucleon loss. The flux of nuclei virith mass number 
A from a source at distance L with initial mass number 
yljni follows simply. 
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Here, la denotes the interaction length of one-nucleon loss 
of the nucleus with mass number a. By employing (II]), 
and summing over all possible mass numbers A (from 1 to 
Aini) different nuclcon loss contribution functions for the 
flux from the source, the total nuclei flux from a single 
source is obtained. 
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this feature, found in Ref. 



and developed further in 



The agreement between the Monte Carlo and analytic 
results, seen in Fig. [T] demonstrates that for negligible 
(<pG) extragalactic fields the local source information 
(spectrum and composition) is essentially encapsulated 
by our simplified analytic model. Indeed, the different 
sets of photo-disintegration cross-sections used for the 
Monte Carlo and analytic descriptions, for which [35] and 
|26j have been used respectively, demonstrate further the 
robustness of these results. Using this analytic descrip- 
tion, sets of models may be quickly scanned through in 
order to find the optimum set of injection index a, expo- 
nential energy cutoff i?max and source composition values 
able to "best-fit" the Auger results above 10^^ eV. 

Through the application of a goodness-of-fit (GOF) 
test, taking into account the systematic uncertainties in 
the measurement of CR energy, (Xmax), and RMS(Xmax) 
(see Appendix^]), whose sizes are shown in Fig. [11 we find 
that, in agreement with our earlier investigations |:16 , a 
intermediate-heavy composition {A > 20), hard spectral 
indices (a < 2) and intermediate type cutoff energies 
(£^Fc,inax ~ 10^"^ eV) are best able to describe the cur- 
rent data. Contour plots showing the "best-fit" regions 
for both silicon and iron only type sources are shown in 
Fig. [2] For lighter nuclei, the contour space was found 
to be considerably diminished, with no such contours ex- 
isting for the proton-only scenario, even at the 99% C.L. 
We also indicate in the plots the position of the GZK 
energy Eqzk in terms of i?Fo,max- The fit to the data 
prefers considerably higher cutoff energies than the GZK 
energy. These plots demonstrate that both the spectra 
and composition information both provide new and dif- 
fering constraints on the source spectral parameters, par- 
ticularly for the case of heavy nuclei (iron) type sources. 
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FIG. 2: Contour plots showing both the spectrum and Xmax 
(ie. both (Xmax) and RMS(Xinax)) goodness-of-fit results for 
silicon (upper panel) and iron (lower panel) source composi- 
tions. The spectral fits to the data have been carried out for 
energies above 10^^ eV, and the QGSJET 11 hadronic model 
has been used as the hadronic model needed to obtain the 
theoretical (Xmax) and RMS(Xniax) values. 



III. THE EFFECT OF THE FINITE DISTANCE 
TO THE NEAREST SOURCE 



As an illustrative example, using typical "best-fit" 
model parameters motivated from our GOF results in 
the previous section, we show explicitly in Fig. |3] the al- 
teration to the flux due to the failure of the homogeneous 
approximation. The lower (upper) panel in Fig. p^ show 
how the arriving flux from an ensemble of sources emit- 
ting iron (silicon) type UHE CR is altered due to a non- 
zero distance to the nearest source. These results demon- 
strate the introduction of a very strong cutoff feature due 
to a non-zero distance to the first source, with this feature 
occurring at lower energy for intermediate nuclei such as 
silicon, than for heavy nuclei such as iron. It should be 
emphasised that the position (in energy) of the cutoff in- 
troduced by the distance to the nearest source is roughly 
independent of the spectral index (a) and cutoff energy 
(i^Fc.max) of the primary spectrum. However, this result 
is only true when the source's maximum energy sits at 
energies much larger than the cutoff energy introduced by 
the local void of sources. We re-iterate here that these re- 
sults rest on the assumption that the composition above 
energies where it has been measured (lO"'^^'^ eV), con- 
tinues to consist of intermediate-heavy nuclei up to the 
highest observed energy flux measurements (10^°'^ eV). 



To put these results on a more general footing, we show 
in Fig. |4] the alteration to the GOF contour plots as the 
distance to the nearest source in increased. We find that 
for both silicon-only and iron-only source scenarios, the 
99% C.L. undergoes a rapid decrease in size for minimum 
source distances in the range 9 — 27 Mpc and 27 — 81 Mpc 
respectively. Furthermore, for larger local voids, the re- 
maining i^Fo.max valucs are larger than those typically 
considered feasible for candidate UHE CR sources. By 
imposing the constraint that the cutoff energy, -BFc.max, 
sits below 10^^ eV, the subsequent upper limit on the 
source distance of 60 Mpc and 80 Mpc are found for sil- 
icon and iron type source respectively, at the 99% C.L. 
This reaffirms the result shown in Fig. |3j that a large 
>81 Mpc local void in the source distribution would lead 
to difficulties in finding self-consistent fits to the spectral, 
(-'i^max), and RMS(Xmax) data of Auger. 



The number density of local sources can not be much 
smaller than 10~^ Mpc~^ as can be estimated, e.g., from 
the absence of "repeaters" in CR data [27l [28j. It is 
hence expected that the fiux of CRs from sources much 
larger than a few 100 Mpc can be well approximated by 
a homogeneous distribution of sources. However, the fi- 
nite (non-zero) distance to the nearest source of UHE CR 
nuclei leads to a breakdown of the homogeneous source 
distribution spectrum result at the highest energies, as 
already indicated in Fig. [ij We here demonstrate the 
effect introduced by this local void of UHE CR sources 
using both concrete examples as well as a statistical ap- 
proach. 



We have assumed throughout this work so far that 
the extragalactic magnetic field strength is negligible. 
Though the overall ffux from the total ensemble of shells 
does not vary when an isotropic and homogeneous distri- 
bution of magnetic fields is introduced, the components 
arriving from different source shells will change. Further- 
more, with an aim to investigate the possibility of a lack 
of nearby UHE CR sources, it is necessary that we take 
into account extragalactic magnetic field effects. In the 
following section we investigate these effects with the aim 
of making our conclusions more general. 
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FIG. 3; The arriving UHE CR flux from a uniform distribu- 
tion of UHE CR sources emitting a pure silicon (top), and 
pure iron (bottom) composition. A source spectral index (q) 
and cutoff energy (-Bpe.max) of 1.6 and 10^°'^ eV respectively 
for silicon and 1.8 and 10'^'^ eV respectively for iron have been 
adopted, as motivated by our goodness-of-fit results in the 
previous section. 
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FIG. 4: Plots showing the disappearance of the 99% GOF 
contours as the distance to the first source is increased up to 
81 Mpc. 



IV. THE EFFECT OF EXTRAGALACTIC 
MAGNETIC FIELDS 



In this section we investigate the effect these fields have 
on the arriving UHE CR flux and composition. In our 
calculations, the extragalactic magnetic field is always 
assumed to have a coherence length of 1 Mpc. Thus, 
within each of the magnetic patches (cells) , the field con- 
tains a uniform component whose direction is assumed 
to be orientated independently (randomly) to that of its 
corresponding orientation in the neighbouring cells. For 
particles with Larmor radii smaller than the magnetic 
patch coherence size {e.g. for iron nuclei with energies 
< 10^^ eV in a nG field), we assume a power law distri- 
bution of magnetic turbulence of the form P{k) ex k~'' 
(where k = 27r/A), for which q = 5/3 corresponds to 
a Kolmogorov-type spectrum, q = 3/2 corresponds to a 
Kraichnan-type spectrum. This angular (and eventually 
spatial) diffusion of the particles is treated following a 



method very similar to that described in [8j, which we 
refer to as the "delta-approximation" method (described 
in further detail in Appendix^]). A Kolmogorov-type de- 
scription of the magnetic field turbulence spectrum is as- 
sumed here. Ultra-high energy CRs in the simulation 
were considered to have arrived once they reached a dis- 
tance of 100 kpc from Earth (this length scale being cho- 
sen to be smaller than both the corresponding particle 
loss length and gyro-radius). 

Provided that a continuous distribution of sources ex- 
ists on all scales, and that the magnetic fields are ho- 
mogeneous and isotropic, extragalactic magnetic fields 
have no overall effect on the arriving fiux (29j . How- 
ever, since a minimum source distance scale must exist, 
a "magnetic horizon" is expected, with the flux from the 
nearby sources being prevented from arriving to us be- 
low a given energy j30j . Furthermore, the contribution 
of sub- "magnetic horizon" source shells can be increased 
by the presence of extragalactic magnetic fields, alter- 
ing somewhat the flux arriving from the different source 
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FIG. 5: A breakdown of the arriving fiux from shells of UHE 
CR Fe emitting sources with an extragalactic magnetic field of 
strength (top) O.I nG and (bottom) I nG. Coherence lengths 
of I Mpc have been assumed for these calculations, which 
were obtained using a Monte Carlo description for UHE CR 
nuclei propagation. 
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shells. 

As seen in Fig. [5] the presence of extragalactic mag- 
netic fields enhances the role played by local regions of 
sources on the arriving UHE CR flux. As envisaged from 
these results, and demonstrated in [31], a low energy cut- 
off of the arriving flux from distant shells of sources is in- 
troduced by the combined effect of a lack of local sources 
and the presence of non- negligible extragalactic magnetic 
fields. Below this cutoff energy, the UHE CR flux from 
even the nearest by sources is suppressed. 

Using the results shown in Fig. [5] the effect introduced 
into the results of section |HI| by the presence of a non- 
negligible (>pG) extragalactic magnetic field may be ad- 
dressed. The dominant effect of these magnetic fields 
in combination with an absence of local sources is their 
alteration of the arriving composition, as demonstrated 
explicitly for 0.1 nG and 1 nG extragalactic magnetic 
fields in Fig. |6] Indeed, with a better handle on both the 
composition and distance to the nearest source, the com- 
position may also provide a valuable probe of the local 
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FIG. 6: The combined effect on the iron spectrum (top panel) 
and composition (middle and bottom panel) due to a local 
void of sources (L> 27 Mpc) and a non-negligible extragalac- 
tic magnetic field (O.I and I nG). For these results, an iron 
only source composition with spectral index a — 1.8) and 



cutoff energy E-Fc,n 



10 eV, have been used. 



V. SUMMARY 



Following the motivation that UHE CRs consist of 
heavy nuclei, whose sources have been suggested to be 
local [321 133] J we have here looked closer at the require- 
ments on their source distribution. In this work we ob- 
tained concrete quantitative constraints on the UHE CR 
source population. Making explicit use of the recently 
provided Auger spectral and shower composition results, 
along with detailed UHE CR nuclei modeling, we inves- 
tigated whether consistency may be found with this data 
using single source composition models. 

For the case of negligible extragalactic magnetic fields, 
we have demonstrated that a simplified analytic descrip- 
tion agrees well with the spectrum and composition re- 
sults obtained from the complete Monte Carlo descrip- 
tion. Utilising this analytic description, we made a scan 
over the source spectral index and exponential energy 
cutoff, for a single source composition scenario, to obtain 
the goodness-of-fit contours. Taking into account the sys- 
tematic errors in the fiux and composition measurements, 
we found that hard (a < 2) source spectral indices and 
intermediate cutoff energies (i?Fc,max ~ 10^°'^ — 10^^ eV) 
for intermediate-to-heavy nuclei could provide a good fit 
to the full set of Auger UHE CR measurements above 
10^^ eV. 

Through the consideration of shells of UHE CR 
sources, we investigated the proximity of the UHE 
CR nuclei sources required if the presently observed 
trend of an increasingly heavy composition continues up 
to the highest energies observed by the ground array 
^2^020.2 gY^^ gy varying the size of the local void up to the 
nearest source, we investigated how detrimental the ef- 
fect of this was on the goodness-of-fits contours. Provided 
that the sources maximum energy lay below 10^^ eV, we 
found that the nearest sources had to be within 60 Mpc 
and 80 Mpc for silicon and iron only sources respectively. 

If, however, extragalactic magnetic fields are suffi- 
ciently strong (>pG), the arriving fiux from the different 
nuclei source shells is altered considerably. Though this 
effect is weakest at the highest energies, a local void of 
sources scenario with such an intervening field, may al- 
ter the arriving fiux at energies below the cutoff feature, 
and thus alter the goodness-of-fit contour landscapes ob- 
tained in Fig.[2J However such (<nG strength) fields are 
found to be unable to alter significantly our upper bound 
on the nearest source distance. 

The requirement for local candidate objects able to 
satisfy the Hillas criterion [53] for > 10^° eV nuclei, 
whilst at the same time not disintegrating these nuclei 
during the acceleration process, places a non-trivial dual 
condition on the source environment. At first consid- 
eration both the required proximity of the sources and 
their required nuclei tolerant radiation fields would seem 
to exclude GRBs as viable candidates. However, recent 
motivations for a Galactic GRB-type origin [55] and the 



possibility that heavy nuclear species are produced dur- 
ing the GRB explosion [36] demonstrate that the case 
for a GRB origin remains viable. Due to the dependence 
of photo-disintegration rates on the radiation field level, 
acceleration sites far away from the central engine are 
naturally favoured. In this regards, acceleration within 
AGN jets or their radio lobes [37] are able to satisfy the 
above mentioned dual condition. 

With regards energetics, the UHE CR source luminos- 
ity density required to power the UHE CR population 
above lO^^ eV is ~ 5 x lO'*'' erg Mpc'^ yr"! [38]. For 
the hard sources spectra motivated in this paper, approx- 
imately the same source luminosity density would exist 
up to the cutoff energy. Thus an average local luminos- 
ity per source of ~ W^^ /Ns erg s~^ is required within 
the local 60 Mpc region, where Ng is the number of con- 
tributing sources. Such a luminosity is roughly l/iVg that 
of the 20-40 keV X-ray luminosity of local AGN [31130], 
of which ~10 sit within 60 Mpc from Earth. Whether it 
is reasonable for these local AGN to channel such a large 
fraction of their non-thermal luminosity into UHE CR 
fiux, and whether a sufficiently enhanced nuclear com- 
position is able to be achieved within their acceleration 
sites, however, remains unclear |41ll42j . 

Our results demonstrate that exciting consequences 
follow from the intermediate-to-heavy nuclei component 
uncovered by Auger measurements. With nuclei photo- 
disintegration inevitably occurring during propagation, 
tough constraints are placed on the source proximity 
and environment. Furthermore, future Auger spectral 
and composition measurements are anticipated to soon 
tighten these constraints. With few candidate sources 
within the present upper bound distance suggested (< 
60-80 Mpc), the puzzle as to the UHE CR origin both 
remains and becomes even more intriguing. 
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Appendix A: Goodness-of-Fit Test 



We perform a goodness of fit (GOF) test of the com- 
patibility of the Auger data with a given model following 
Ref . ^43] . For a fixed source composition of nuclei we vary 
the universal spectral index a of the emission spectrum 
and the maximal rigidity cutoff -Epcmax; described in [45] 
and [46]. 



Given the acceptance A„ (in units of area per unit time 
per unit solid angle) of the experiment in bin n centered 
at energy En with bin width A„, the number of expected 
events is 



^„ (a, £'Fc,max, A/", 5i 



E+ 



AnAfJdEY^^iE), (Al) 



£,7 



where dNi /dE is the flux of nuclei Z arriving at the de- 
tector and the boundaries are E^ = £'„(! + 6e) ± A„/2. 
The parameter Se in Eq. ( Al ) is a fractional energy- 



scale shift that takes into account the uncertainty in the 
energy-scale and Af is the normalization of the source 
luminosity. 

A fraction e„ of the expected /i„ events per bin passes 
the quality cuts for the distribution of shower maxima 
X. This distribution is only known by its first two mo- 
ments, the average shower maximum {X)n and its root- 
mean-square (RMS) AXn [2- For the GOF we hence 
divide the distribution of shower maxima X into three 
intervals (X~„,X+„), such that each interval contains 
1/3 of the total number of observed events. For a Gaus- 
sian distribution this corresponds to a central bin of size 
Xo±„ ~ {X)n{l + S^x)) ± 0.43AX„(1 + Snus)) and two 
bins containing all other events left (m = — 1) and right 
(m = 1) to it. Here, we introduce the fractional un- 
certainties 6(^x) and iJrms of the mean and RMS of the 
shower maximum, respectively. 

We emphasise that the binning method of the shower 
maxima that we employ can only as an approximation. 
As a check of this approach we try to reconstruct {X) 
and AX from this X-binning by a simple x^-fit. The 
number of events passing the quality cuts is known from 
Fig. 2 of Ref . [9, . Figure It] shows the reconstructed la 
contours in comparison to the statistical uncertainty of 
the measurement extracted from Fig. 2 and 3 of Ref. [Sj . 
The contours are consistent with the data and reproduce 
the right order of magnitude of the statistical uncertainty. 

The predicted distribution of shower maxima can be 
determined by hadronic interactions models such as 
QGSJET 11 [lO]. We assume in the following that the 
distribution in the n-th energy bin can be approximated 
by a Gaussian distribution with mean X„ and RMS 
cr„. The expected number of events in the X sub-bins 
m = —1, 0, 1 is hence of the form 

^J'm.n{ol■^ £^Fc,max, A/", 5ei & {X) : (^RMs) 



e„/i„(a, E^Fc.maxjA/", (5b) / dX 



2^ 



27rCTr, 



(A2) 



X' 



The probability P„(A^tot,n) of observing iVtot,n events 
in the n-th energy bin follows a Poisson distribution 



f{Ntot,n,IJ'n) with mean /x„. The e„7Vtot,n events pass- 
ing the quality cuts distribute between the three X-bins 
following a Poisson distribution with mean fim.n subject 
to the constraint enNtot.n = A^-i,n + A^o.n + A^i.n- The 
conditional probability is 



P;(iV_i,„,iVo,„,^i,„) 



(A3) 



Hence the total probability of observing a set of events 
{N} is given as 

P{{N}; a, i^Fcanax, A/", (5(x) , Sb.ms,Se) 

- n ^"(^tot,„) n ^^^-1,^' ^0,^, Ni,i) ■ (A4) 

n i 

Here the first product runs over all energy bins used for 
the CR spectrum and the second product runs over all 
bins where additional information about the elongation 
rate distribution is available. 

As a next step we marginalize over the experimental 
systematic uncertainty in the energy scale Se, shower 



maximum 



5(x) and root-mean-square iJrms as well as 
the normalization A/". We assume a flat prior with 
\Se\ < 20%, |(5(x)l < 2% and |5r,msI < 10%. Marginal- 
ization is done by maximizing the probability of the ac- 
tual experimental observation. We hence define 

P({7V};a,^Fc,max) 

= P({Af};a, £'Fe,max,A/', (5(x),(5rms,^£;), (A5) 
where the set of parameters {J^,S^x),SrmS7Se} maxi- 



mizes the probability ( A5 1 for the experimental result 
jjYcxpj.^ consisting of the total events per bin N^^^^. as 
well as N^^^ = enN^o^^^/'^t if available, which follows 
from the construction of the X-binning. 

The transition between the galactic and extragalactic 
component of UHE CRs is uncertain. A natural can- 
didate is the CR "ankle" that is observed in the Auger 
data at about 10^^'^ eV. We assume that the "contami- 
nation" by a low-energy galactic component has become 
negligible beyond 10^^ eV and include only Auger data 
above this threshold for our GOF. For consistency we 
have to check that the extragalactic component does not 
over-shoot the CR data at lower energies. Hence, in the 
marginalization procedure we also impose a prior on the 
normalization Af by requiring that the model spectra do 
not exceed the Auger data below the first bin used in the 
fit by more than three standard deviations. 

Finally, the model (a, -BFo,max) is compatible with the 
experimental results at a given GOF if 

J2 P({A^};a,^Fe.max)<GOF. (A6) 

P({Af})>P({JV<=^P}) 
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FIG. 7: The statistical uncertainty of the measured mean 
and RMS of the shower maximum from Ref. 9 compared to 
the uncertainty of the reconstruction from the three X-bins 
introduced for the goodness of fit. 



Technically, this calculation is performed by generating 
a large number {iV'^P} of replica experiments following 
the probability distribution ( A5 1 and by imposing that 
a fraction F with P({7V''=p}) > P({7V'='^p}) satisfies F < 
GOF. 



We prepare a turbulent magnetic field region, contain- 
ing an ensemble of waves with wavelengths between Amin 
and A,„ax- It is assumed that a guiding mean- field, Bq, 
lies along the z-direction, and that the sum of the energy 
within the turbulent field is equal to the total energy 
within the guiding field {Bq — J2n^^n)- We describe 
the turbulence using 50 isotropic waves, with 10 waves 
per decade (i.e. Amax/Amin = 10^). The Bulirsch-Stoer 
method was used to track the particles in the field. 

A simplified description of turbulent propagation, how- 
ever, may also be obtained through a scattering descrip- 
tion of UHE CR interaction with the turbulent field. In 
this description, UHE CR stochastically interact with 
scattering centers whose density is dictated by the energy 
density in the turbulence power spectrum at the wave- 
length matching the gyro-radius of the UHE CR. We re- 
fer to this description as the "delta-approximation" case, 
since it assumes only resonant scattering. The resonant 
scattering angle of this process is given by 



Ae 



i?L > Lcoh , 



(B2) 



where i?L is the particle's Larmor radius and Lcoh is 
the magnetic field's coherence length. The corresponding 
scattering length is then given by 



i?scatt(A6l)=aL< 



J(i) 



2-9 



1 



i?L < icoh 
i?L > Lcoh 



(B3) 



Appendix B: A Simplified Description for UHE CR 
Propagation Through Turbulent Fields 



We here describe the results of an investigation we car- 
ried out comparing two different descriptions for the dif- 
fusive propagation process of UHE CR in turbulent mag- 
netic fields. 

As discussed in [U] , a turbulent field may be generated 
by adding plane waves whose amplitudes are dictated by 
the magnetic power spectrum. 



B 



Y^SBr^ 



(Bl) 



where (5B„ = ^n^ne'-*''"'^'^'^"), with x' = xcosa - 
ysina, a describing the orientation of the wave (ran- 
domly chosen), £„ describing the (randomly chosen) po- 
larisation of the wave (ie. over each wavelength ^„ de- 
scribes an ellipse in the y' — z' plane), A^ describing 
the amplitude of the wave, fc„ (= 27r/A„) describing the 
wavenumber, and 0„ is a random phase. The amplitudes 
of the different waves are given a power law distribution 
of the form An (x kn'' ■ For the purpose of our compar- 
ison we assume a Kolmogorov spectrum (i.e. q — 5/3). 



where the factor a describes how far this scattering is 
from the Bohm regime (i.e. a = 1 at the Bohm limit, 
giving rise to approximately one scattering per gyro- 
radius) for particles whose gyro-radii are at the scale 
Amax- The factor (2 — q) may be thought to describe 
the exponent in the energy dependence of the density of 
(resonant turbulence) scatterers. In order to scale down 
this result to allow us to probe the low angle scattering 
regime, we employ a small angle scattering length of size 
Se, Rscttm = i5e/AeyRsc,tt{Ae). Note, the intrinsic 
scattering in angle introduced through each photodisinte- 
gration process is always heavily sub-dominant, and can 
be safely ignored. 

We compare the growth of spread of particles in the 
mean-field direction (Az) as a function of time for both 
the simplified description and that of the full turbulent 
field description in Fig. [8] for the case of a = 1. For this 
simulation, the particles were all started from the same 
position in the static turbulent field with an isotropic 
distribution of directions. 

With good agreement found between the "full power 
spectrum" and the "delta-approximation" descriptions 
we conclude that the "delta-approximation" method is 
able to provide an accurate description for UHE CR prop- 
agation in extragalactic magnetic fields. 
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FIG. 8: A plot comparing the spreading of particles along 
the mean-field direction during propagation for both the 
full turbulent field description and the simplified "delta- 
approximation". These results were obtained for 10^* eV 
energy UHE CR iron nuclei in 0.1 uG, 1 nG, and 10 nG ex- 
tragalactic magnetic field with a 1 Mpc coherence length. 
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